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We have developed a subtractive renormalization method with which we can evaluate nucleon- 
nucleon (NN) scattering phase shifts produced by the NN potential obtained at leading, next- 
to-leading, and next-to-next-to-leading order (NNLO) in chiral effective theory (#ET). In this 
method the low-energy constants associated with short-distance NN physics are eliminated from 
the Lippmann-Schwinger equation (LSE) for the NN t-matrix, in favor of physical observables. 
This allows us to straightforwardly compute scattering phase shifts for ultra-violet cutoffs of at 
least 10 GeV. We then perform detailed analyses of the maximum cutoff at which the use of a 
%ET NN potential in the LSE makes sense. 
Specifically, we show that: 

(a) our subtractive renormalization technique reproduces known results for the LO potential, in 
both S- and P- waves; 

(b) a parameterization of short-distance physics in the NNLO potential in terms of an energy- 
dependent contact term creates scattering resonances and shallow bound states in S-wave channels 
once cutoffs larger than 1 GeV are considered; 

(c) the more conventional momentum-dependent contact term in the NNLO potential has prob- 
lems of its own at cutoffs larger than 1 GeV; 

(d) the NNLO potential yields P-wave phase shifts that have significant dependence on renormal- 
ization point. 

(e) for cutoffs smaller than 1 GeV, using spectral-function regularization for the long-distance part 
of the potential produces results that vary with the cutoff and depend on the renormalization point 
less than if dimensional regularization is employed to compute the two-pion-exchange graphs. 
Based on all these results we conclude that, once cutoffs larger than the chiral-symmetry breaking 
scale are employed, iteration of the two-pion-exchange piece of the JET NN potential in the LSE 
does not satisfy all of the criteria required for successful renormalization of the problem. 
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1. Introduction 

Chiral perturbation theory (#PT) is an effective field theory (EFT) that enables calculation 
in the non-perturbative region of QCD. The use of EFT in nuclear systems conveys two major 
advantages. First, we gain physical insight into the behavior of the strong interaction at large 
distances, which is difficult to calculate ab initio from QCD. Second, our calculations become more 
accurate order by order, thus allowing systematic error control. An EFT becomes most powerful 
when a clear and large seperation between the low- and high-energy scale in the problem is possible. 
It has been shown that #PT is quite successful in the low-energy (< 1 GeV) mesonic sector. In 
principle, the same theory should work in the low-energy nucleon-nucleon (NN) sector as well. 

However, it has been almost two decades since #PT was first applied to the problem of NN 
system||l|], and difficulties still remain. Standard #PT power counting, which would predict no 
bound state for the deuteron, does not apply to the NN problem, because of infrared enhancements 
of the NN interaction. A nonperturbative treatment of at least part of the NN interaction is thus a 
necessary ingredient. One needs to either iterate the NN potential computed from #PT using the 
Lippmann-Schwinger equation (LSE) || |, |, |, g g g |fi| |ll| , [T|, |l3[ 0, [T|, [l| 0, |E]], 



or determine which part of the potential can be treated in perturbation theory Q2Q , 22J , |22|, |23, |24j, 



25| , |26| |27p. So far there is no consensus as to which of these two alternatives is superior. 

Here, we adopt the former approach. In this "chiral effective theory" (#ET) the behavior of the 
#PT potential at high momentum necessitates that a cutoff A be placed on the momenta in the LSE. 
It is then natural to ask what values of A can be used, if renormalization is to be successfully carried 
out. Difficulties in answering this question have both a technical part, i.e. it is hard to perform a 
fit for all unknown low energy constants at high cutoffs because of "fine-tuning"; and conceptual 
problems regarding what a successful renormalization is, e.g. is it sufficient that observables be 
(approximately) cutoff-independent? (See, e.g., Ref . J28|] , for a recent discussion.) 

In Sec. ^ we outline a subtractive renormalization technique that solves the "fine-tuning" prob- 
lem. This technique thus allows us to assess how well %ET at large cutoffs satisfies criteria for 
successful renormalization. We have used subtractive renormalization to calculate the NN scatter- 
ing amplitude obtained by using leading-order (LO), next-to-leading-order (NLO), and NNLO #PT 
NN potentials in the LSE, for cutoffs up to A = 19 GeV. We show some results of these calculations 
in Sec. [3| and examine the conditions under which %ET is really improved, order by order, after 
renormalization. We do this for both the dimensionally-regularized (DR) and spectral-function- 



regularized (SFR) [29] #PT potentials, and consider both energy and momentum-dependent con- 



tact terms. More details regarding all these methods and results can be found in Refs. 31), 32], 

2. Main ideas of subtractive renormalization 

The main idea of our subtraction method is to construct the fully off-shell partial-wave ^-matrix 
from the knowledge of the long-range part of the potential and the on-shell value of the ^-matrix 
for zero energy [30, 31, 32, 33, 0]. The partial- wave LSE is given by 



t n {p ,p;E)=v n (p ,P) + y-M T —. . (2.1) 

yf 71 Jq p^ + ie- p" L 
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Where Pq/M = E is the center-of-momentum (cm.) energy and A the cutoff parameter. The 
incoming (outgoing) angular momenta are indicated by / (/')■ The potentials are defined as: 



v in (p',p) = vf,f( P ',p) +Q, l p n 'p l f(p,p , ;E 



(2.2) 



where p(p') indicates the incoming (outgoing) momentum in the cm. frame, vffl is the long-range 
potential that is operative in this channel. Cmp' 1 ' p l f represents the contact interaction, where / can 
be energy or momentum-dependent. First we consider the case / = 1, which is the simplest contact 
term in a given partial wave. To relate the ?-matrix to a physical quantity, a generalized scattering 
length for arbitrary angular momenta / and V can be defined as [35] ^ = liirn^o ''''^f^ , where 
for /' = I = the usual definition, ^ = ? o(0,0;0), i s obtained. Dividing the partial-wave LSE, 



Eq. (|2. 1[), by p n p l we obtain 



tfUp'^E) _v s / l (p',p) + y2 



iV I 

p" p l 



p'l'p' 



M 



jf 7C p"'p l 



*dp" p" 2 vf l i ll {p',p")tf l f [ {p",p;E 
pl + ie - p" 2 



(2.3) 



Since vf,f(p',p) ~ p n p l , Eq. (23) is general and can be applied to any partial wave. 

In the following we concentrate on P-waves (/ = /'= 1). (The corresponding argument for S- 
waves, in the case that we have the standard LO contact interaction of #PT with I = I' = 0, /= 1, is 
analogous, but more straightforward, as division by a factor of p'k is not necessary there.) Consider 
the half-shell and on-shell ^-matrices at E = 0: 
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Subtracting Eq. (| 

~tvi(p',k;0) 



lim 

it->o 



from Eq. (2A) and multiplying both sides by p' cancels the unknown Cyf. 
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Here we have used that for P-waves lim^o 
The only unknown in Eq. (E1J) is lim^o 



fu(U;Q) 
kk 



j^-. The above limits are well-defined. 



k 



, which can be solved by standard techniques. 



The next step is to apply the same idea again to obtain 



t,'i(p,p';Q) 



and hence, ?/'/(/?, p';0). We 



then proceed to calculate the on-shell ?-matrix and the phase shifts using resolvent identities that 
connect the operator t(E) to the operator f(0). Those details are laid out in Refs. [34, 30, 39]. 

Next, we consider an energy-dependent contact term in S-waves. We take /Coo = A + jE. This 
is the contact term up to NLO and NNLO in #ET for the 1 Sq channel. To simplify the presentation, 
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we adopt the following operator notation for the LSE 

t(E)=X + yE + v LR +[X + yE + v LR ] g (E)t(E), (2.7) 

where go(E) is the free resolvent of the LSE. Setting E = in Eq. (2.7) leads to t(0) = X + vlr + 
[k + vlr] go(0) f (0), which contains only one unknown, A. Therefore, the matrix element t(p',p;0) 
can be obtained from one experimental datum, here the NN scattering length oiq. After applying the 
same idea to obtain t(E*) from the phase shifts at an arbitrary energy, and using properties of the 
LSE to eliminate y, we have 

t(E) +t(0)[ go (0)-g (E)]t(E) + §;{t(0) - [1 +t(0)g (0)] oct(E*)} g (E)t(E) 

= ~ ^) t{0) + ^°(°)] at( - E *^ (2 - 8) 



where a = [1 +t(E*)g (E*)}- 1 . With t(0) and t(E*) known, Eq. (|ZJ|) is an equation for t(E) and 



can be solved by standard methods with detail given in Ref. []32[]. 
Finally, for (more complicated) contact terms such as 

we can solve the problem by combining the above methods, i.e., use the first subtraction to elim- 
inate A, and then relate t(E*) to t(E). For coupled channels, we can apply the idea of dividing 
p' 1 ' p l in LSE to eliminate X t p 2 . However, in the case of the momentum-dependent S-wave contact 
terms, we need to perform one fitting to eliminate the unknown constant C<i. The inputs needed are: 

for case (A); oq, and an additional data to perform the fitting; 

for case (B); ao, (X20, and an additional data to perform the fitting; and 

for case (C); ao, OI20 and phase shift at an arbitrary energy 8(E*). 1 



3. Results and discussion 



In this section we present our results in P-waves and S-waves to demonstrate the following: 

(1) . Our subtractive renormalization scheme generates results equivalent to the conventional 
"fitting" method, with a direct input of physical observables. 

(2) . The energy-dependent contact term produces phase shifts that oscillate with respect to A. 

(3) . Whether a contact term (or contact terms) is needed for a cutoff independent result is 
exactly determined by the (coordinate-space) singularity structure of the potential as r — > 0. 

(4) . Cutoff independence in the phase shift does not neccessarily mean the results are renormalization- 
point independent. Both properties are necessary conditions for successful renormalization. 

(5) . In general, there is a highest cutoff A c w 1 (2) GeV in the LSE one can adopt for the 
NNLO DR (SFR) TPE, before the results start to exhibit problems. 

To show point (1), we compare the LO S-wave phase shift obtained from the conventional 
"fitting" method to our subtractive scheme in the left-hand side of Fig. |l| Here the potential is 
the one-pion-exchange (OPE) plus a constant contact term. Fig. [l] shows that the results obtained 

'The only restriction on E* is that it must be within the domain of validity of our theory. 
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Figure 1: (Color online) The comparison of two renormalization methods for the lowest NN 1 So & 3 Si- 3 Di 
phase shifts (left two panels) and the bound-state wavefunctions (right two panels). y/b(p) (u(r)) is the 3 S] 
wavefunction and ^(p) ( w ( r )) denotes the 3 Di wave in momentum (coordinate) space. Here A = 50 GeV 
is used. The dotted lines indicate the corresponding results obtained with the CD-Bonn potential 

by these two methods agree with each other within a relative difference of 2%. (Due to numerical 
effects this difference is amplified in E\. E\ is well known for its sensitivity of the value of the 
unknown constant in the contact term, and we only adjusted ao up to certain precision when per- 
forming the "fitting" method.) We have verified that the off-shell t-matrices we obtain also agree 
with those obtained from the conventional "fitting" method to the same accuracy. Therefore, our 
subtraction method is as valid as the conventional fitting method. Our method can be applied to 
bound-state calculations too. The right-hand side of Fig. [T] shows deuteron wavefunctions obtained 
from our subtractive method, which are quite close to those obtained from the CD-Bonn potential. 



> 




2000 4000 6000 8000 1000012000140001600018000 
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Figure 2: The 1 S NN phase shift at T, uh = 10 
(upper panel) and 100 (lower panel) MeV as a 
function of A. The results are obtained using 
the DR NNLO TPE with an energy-dependent 
contact term via our subtractive renormaliza- 
tion. 



For point (2), we associate the DR NNLO TPE with the energy-dependent contact term and 
plot the 1 So phase shifts versus A in Fig. 2. The phase shifts show an oscillatory behvaior as a 
function of A. A similar oscillatory feature is observed in the 3 Si — 3 D] channel. (We use ao, 



5 



Subtractive renormalization of the NN system 



C.-J. Yang 





'Pi 


3 Po 


3 Pi 


3 P 2 


Table I: Singularity structure of the 


OPE 


U 


R 


u 


R 


long-range potentials v\f. Here "U" 


NLO (DR) 


U 


U 


R 


R 


("R") means that Vjf is repulsive (at- 


NNLO (DR) 


R 


R 


R 


* 


tractive) at r — » 0. The * indicates that 


NLO (DR) + NNLO (SFR) 


U 


U 


R 


R 


both eigenpotentials in the 3 P 2 — 3 F2 


NNLO (SFR) 


U 


U 


R 


R 


channel are attractive. 



5 [deg] 



20 
15 
10 
5 

-5 
-10 



Unrenormalized 

8 [deg] 



Re normalized 



- Ill 1 
II 

- II 
II 
II 
II 

T '1 


1 1 1 

'p, 
t 


— SI 


_ r 


• ■* tv- ■ ■ 


■f- 






— 1-^ 




1000 2000 3000 4000 
A [MeV] 




1000 2000 3000 4000 
A [MeV] 



1000 2000 
A [MeV] 



1 000 2000 
A [MeV] 



Figure 3: (Color online) The un-renormalized v.s. renormalized NN P-wave phase shifts at 7) a b =10 MeV 
(un-renormalized) and 100 MeV (renormalized) as a function of A for various ^PT potentials: DR NLO, 
black dotted line; DR NNLO, red dashed line; SFR NNLO, solid green line. For the renormalized case, 
the input a\\ were adjusted at each cutoff to give the best fit to the Nijmegen analysis [^6|| in the region 
T lab < 100 MeV. 



8(E*) and a 2 o (for the triplet) as the input to generate the results with E* = 1.4 (10) MeV for the 
singlet (triplet) channel.) This phenomenon is caused by the resonance state created by the energy- 
dependent potential. We emphasise that the first place where the phase shifts diverge is at cutoff 
A« 1 (1.2) GeV for the singlet (triplet) channel. 

(3) involves the short distance (r — > 0) behavior of vf,f in the coordinate space. We calculate 
the r — > behavior analytically for various P-waves potentials and list them in Table I. At the same 
time, we plot the un-renormalized v.s. renormalized phase shifts at Tig}, = 10 (100) MeV for TPE 
up to DR NLO, DR NNLO and the SFR TPE up to NNLO in Fig. 3. (Here, and throughout 
this paper, we adopt an intrinsic cutoff A = 800 MeV for the SFR TPE.) For the renormalized 
cases, the contact term has the form Cf,jp' 1 p l . Comparing the un-renormalized v.s. renormalized 
case indicates whether a contact term is needed for the phase shift to be stable with respect to A. 
This is exactly determined by the r — > structure listed in Table I. If the potential is singular and 
attractive for r — > (denoted as"R" in Table I), then the contact term is required. If it is not ("U" 
in Table I) then the phase shifts will have a stable A — > 00 limit even in the absence of a contact 
term (see also Refs. [0, [Tj§ 0). The 3 P 2 - 3 F 2 channel for DR NNLO is a special case, since 
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the coupled-channels potential has two attractive singular eigenpotentials in the r — > limit, and so 
one subtraction is not sufficient to make phase shifts independent of A in this channel. 
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Figure 4: (Color online) The % NN 
phase shift as a function of the lab. ki- 
netic energy for various A. The results 
are obtained with the DR (left two pan- 
els) or SFR TPE (right two panels) up to 
NNLO with momentum-dependent con- 
tact terms. We use oq = —23.7 fm as in- 
put and then perform a fit to either the 
effective range ro = 2.7 fm (solid black 
line) or the phase shift at T[ a t, = 200 MeV 
(dashed red line). The phase shifts 
are denoted by open triangles. 



(4) To see why a cutoff-independent result in the phase shift is not neccessarily renormalization- 
point-independent, we plot the 1 Sq phase shift obtained with the DR and SFR NNLO TPE along 
with the momentum-dependent contact term (denoted as case (A) at the end of the previous section) 
in Fig. 4. As mentioned before, in this case we perform the renormalization by the one-subtraction- 
plus-one-fitting procedure. The results obtained by fitting to the effective range r a or to the phase 
shift at Tiab = 200 MeV are shown. One can see that the two different fit procedures generate 
different results for the same A. This is especially visible at A = 500 and 1000 MeV for the DR 
NNLO TPE, where a resonance-like behavior is present in the latter case when C2 is fitted to ro. 
For values of A not close to these problematic cutoffs the phase shift is almost independent of the 
renormalization point. In contrast, for the SFR TPE, the two different fitting procedures lead to 
almost the same phase shift for A between 700 — 1800 MeV. By switching to the SFR TPE, we 
achieve renormalization-point-independence for a wider range of A. 



DR NNLO 



A=600 
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■ ■ A=800 

- A=900 

— A=1000 
Nij93 




Figure 5: (Color online) The best fit for 
the NN 3 Si— 3 Di phase shifts as a function 
of the laboratory kinetic energy for different 
cutoffs A ranging from 0.6 to 1 GeV. The po- 
tentials employed are the SFR NNLO with 
a momentum-dependent central part of the 
contact term. The values of the Nijmegen 



phase-shifts 1 36 1 are indicated by the open 
triangles. 
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Finally, for point (5), we plot the 3 Si — 3 Di phase shifts in Fig. 5. These results are obtained 
by the SFR or DR TPE up to NNLO plus the momentum-dependent contact term (labeled as (C) 
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in the previous section). For the DR TPE case, the best overall fit already diverges away from 
the Nijmegen analysis in the mixing angle at A = 1000 MeV. This implies that there is a critical 
cutoff A c ~ 1 GeV. Above that we cannot iterate DR NNLO TPE in LSE and obtaine a good 
fit in the 3 Si — 3 Di channel. Moreover, as shown in Fig. 4, at this A c the renormalization-point- 
independence also breaks down for the DR NNLO TPE in the ! So channel. Therefore, we conclude 
that for the DR TPE in S-waves, the highest cutoff one can adopt in the LSE is A c ~ 1 GeV. For 
the P- waves, a detail analysis of the renormalization-point-dependence suggests that A c ~ 1 — 1 .2 



GeV for the DR NNLO TPE[31]. As with the S-waves, if the SFR TPE is adopted, then A c can be 



extented to 2 GeV before similar problems appear. 

4. Summary and Conclusions 

We developed a subtractive renormalization scheme for #ET NN potentials which allows us to 
go to an arbitrarily high cutoff in the LSE. Our calculations show that the energy-dependent contact 
term creates scattering resonances and shallow bound states in S-wave channels once cutoffs larger 
than 1 GeV are considered. Momentum-dependent contact terms in the NNLO potential also has 
problems at these cutoffs. We also investigate the singularity structure of the potential and find that 
the LO conclusion presented in Ref . [|l5|] holds up to NNLO. Our analysis in S-waves and P-waves 
shows that the two-pion-exchange potential should not be inserted in the Lippmann-Schwinger 
equation and treated non-perturbatively if cutoffs larger than 1 GeV are employed. 
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